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Delayed pulsar kicks from the emission of sterile neutrinos 
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The observed velocities of pulsars suggest the possibility that sterile neutrinos with mass of several 
keV are emitted from a cooling neutron star. The same sterile neutrinos could constitute all or part 
of cosmological dark matter. The neutrino-driven kicks can exhibit delays depending on the mass 
and the mixing angle, which can be compared with the pulsar data. We discuss the allowed ranges 
of sterile neutrino parameters, consistent with the latest cosmological and X-ray bounds, which can 
explain the pulsar kicks for different delay times. 

PACS numbers: PACS numbers: 97.60.Gb, 14.60.Pq, 97.60.Bw 



Observed velocities of pulsars pj] can be explained by 
an anisotropic emission of sterile neutrinos [1, S 0, Q 
or other light particles @. Sterile neutrinos are firmly 
rooted in particle physics 0], because the gauge singlet 
(right-handed) neutrinos are needed to account for the 
observed neutrino masses in what is called the seesaw 
Lagrangian Q . If some of the gauge singlets turn out to 
be light, they can appear below the electroweak scale as 
sterile neutrinos. There are further hints in favor of this 
intriguing possibility: in addition to explaining the pulsar 
kicks, sterile neutrinos with the same parameters could 
make up the cosmological dark matter [l^, |Tl|, [l^, 
and could play an important role in the formation of 
the first stars This form of "warm" dark mat- 

ter may be in good agreements with observational in- 
ferences regarding the small-scale structure [l^ [l^ . In 
contrast, the active neutrino oscillations cannot explain 
the pulsar kicks, unless they have very large magnetic 
moments 117|, or the mass difference is large enough to 
allow the Mikheev-Smirnov-Wolfenstein ^8*1 resonance at 
density 10" - lO^^ g/cm^ [ll,!!^, which is excluded by 
the current data on neutrino masses. 

The most promising way to discover sterile neutrinos 
is by observing X-ray photons from decays of the relic 
sterile neutrinos, which have lifetimes longer than the 
age of the universe but which can, nevertheless, produce 
a detectable signal [2l[. The X-rays produced by the 
dark matter decays, and the production rate of sterile 
neutrinos in a supernova are both governed by two pa- 
rameters: the sterile neutrino mass and the mixing an- 
gle. It is, therefore, important to understand the allowed 
ranges of these parameters. The range implied by the 
pulsar kicks can help focus the X-ray searches. In Ref. 3, 
the allowed range was discussed for both the resonance 
and the off-resonance production. In this paper we will 
reconsider this range of parameters, apply the present 
constraints from the X-ray observations and the Lyman- 
alpha bounds, and we will also relate the range of sterile 
neutrino masses and mixing angles to the delays in the 
onset of the kick. This delay can have some observable 
consequences and can be determined from the studies of 



the pulsar populations |23 |. 

In applying the cosmological constraints, we distin- 
guish between two different issues: the particle's exis- 
tence and its ability to account for all of dark matter. 
Assuming the standard cosmology (rather than, e.g., the 
low-reheat cosmology p3j). one expects the neutrino os- 
cillations to generate some out-of-equilibrium population 
of dark-matter particles that depends only on the mass 
and the mixing parameters. This population may con- 
stitute only a fraction of dark matter if the mixing angle 
is small enough. However, there are other ways in which 
the relic sterile neutrinos could be produced: for example, 
they can be produced from the inflaton decay or from 
the decay of the Higgs boson at the electroweak scale [ll| . 
The latter scenario produces a population of dark mat- 
ter particles that is considerably "colder" [11] than the 
warm dark matter originating from neutrino oscillations, 
and the amount of dark matter is completely indepen- 
dent from the mixing angle. Nevertheless, regardless of 
any additional production mechanisms, the production 
by oscillations cannot be "turned off" (except for non- 
standard cosmological scenarios 0]). Therefore, there 
exists a robust cosmological bound on the mass and mix- 
ing angle, which is based on the effects of sterile neutrinos 
produced by oscillations, even if they do not make up all 
the dark matter. 

We, therefore, show two exclusion regions in Fig. 1. 
In the solid "excluded region" , the existence of a sterile 
neutrino conflicts with the assumptions of standard cos- 
mology. Below this region, but above the dashed line, the 
particle may exist, but it may not account for all the dark 
matter because of the existsing X-ray bounds [2l|. Fi- 
nally, below the dashed line, the particle may exist and 
may account for all dark matter. As one can see from 
the figure, the resonant mechanism is inconsistent with 
all dark matter being sterile neutrinos, but is a viable 
explanation for the pulsar kicks, as long as the sterile 
neutrinos make up only a part of dark matter. In con- 
trast, the off-resonant mechanism is consistent with all 
of dark matter being in the form of sterile neutrinos. 

In a supernova, the sterile neutrinos are produced in 
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two ways: the active neutrinos can oscillate into the ster- 
ile neutrinos on resonance 0, or ofF-resonance In 
both cases, they escape anisotropically because the elec- 
trons and other fermions in the newly formed neutron 
star are polarized in the magnetic field. Of course, the 
ordinary neutrinos are produced with some anisotropy 
as well, but their production asymmetry is completely 
washed out by the numerous scatterings the neutrinos 
undergo on their way out of the star as they diffuse 
in approximate thermal equilibrium [20| . In contrast, 
the sterile neutrinos escape without scatterings, with the 
emission asymmetry equal their production asymmetry, 
because their scattering cross section is suppressed by the 
small mixing angle. 

In the case of the production off resonance, the allowed 
range of parameters has a direct connection with the time 
delay from the supernova collapse until the onset of the 
kick. In the case of production on resonance, the connec- 
tion is less obvious. While there are many uncertainties 
in the supernova parameters, we think it is of interest 
to show the parametric dependence of the allowed kick 
parameters on time delay. 

Let us briefly summarize the results of Ref. The 
off-resonance production rate of sterile neutrinos is de- 
termined by the mixing angle in matter 9m, which, in 
general is not the same as the mixing angle 9 in vacuum: 



sin^ 29m. = 



(AmV2p)2sin2 261 



( Ato2 /2p)2 sin^ 29 + (Am^ /2p cos 29 ~ V^f ' 

(1) 

where the matter potential Vm is positive (negative) for 
v(y)-, respectively; p is the momentum. For the case of 
oscillations into sterile neutrinos. 
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(SY, - 1 + AY,^ + 2Y,^ + 2Y,J. (2) 



In a core collapse supernova, the initial value of this 
matter potential is Vm — (—0.2... -I- 0.5)Vo, where 
Vq depends on the density p; Vq = Q,p/V2m„ ~ 
3.8 eV(p/10i4gcm-3). 

It was pointed out in Ref. that, in the presence 
of sterile neutrinos, rapid conversions between different 
neutrino flavors can drive the effective potential to its 
stable equilibrium fixed point 

Vm ^ 0. (3) 

This equilibration takes place on a time scale 
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Once the equilibrium is achieved at Vm « 0, the effec- 
tive mixing angle in matter is close to that in vacuum, 
and from this point on the emission of sterile neutrinos 
proceeds at a much higher rate. 
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FIG. 1: The allowed regions for delayed kicks with delays 
from 1 through 5 seconds (assuming the other parameters are 
fixed) are shown by black solid lines marked by the numbers 
representing the delay time in seconds. The exclusion region 
is based on the combination of the X-ray bounds and the 
Lyman-a bounds, according to Palazzo et al. P^]. Here we 
distinguish between the two possibilities: (i) sterile neutrino 
with a given mass and mixing angle may exist, and (ii) sterile 
neutrino with a given mass and mixing angle may constitute 
the entire cosmological dark matter. The former possibility 
is viable for all the points below the solid "excluded region" , 
while the latter is limited from above by the dashed line (see 
discussion in the text). These bounds can be evaded in low- 
reheat cosmologies [231 ]. 



There is a considerable uncertainty in the equilibra- 
tion time given by eq. (O because several parameters are 
functions of time and position in the star and are not 
known precisely. The equilibration does not have to oc- 
cur simultaneously in the entire star. However, although 
the emission of both active and sterile neutrinos is sub- 
ject to much uncertainty, the total energy must be close 
to the initial gravitational energy of the collapsing core, 
3 X 10^"^ erg. It is also known that at least 30% of this 
energy must be carried out of the supernova by ordinary 
neutrinos to explain the observed neutrino signal from 
SN1987A. Finally, if the sterile neutrinos are to explain 
the pulsar kicks, they must carry a non-negligible fraction 
£s of the total energy £tot : 



^^tot 



0.25. 



(5) 



This results in a 1% anisotropy of the overall neutrino 
emission and give the pulsar a kick consistent with ob- 
servations. 

The relative rates of the active and sterile neutrino 
production depend on the mixing angle 9 and the tem- 
peratures in the core (where the sterile neutrinos are 
produced) and the neutrinosphere (from where the ac- 
tive neutrinos are emitted). In addition, the time inter- 
val over which the sterile neutrinos are emitted can be 
shortened if the equilibration of Vm — > takes some non- 
negligible time. Based on the discussion of Ref. one 
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can estimate the ratio (fg/^^tot) = £s/{£s + £v) as follows: 
(ll)^sinH (-^^X ('—^) f^Uo,., (6) 

J — sphere J \ ^ J 

where < 1 is the fraction of enclosed mass of the 
core from which the emission of sterile neutrinos is effi- 
cient, and /d.o.f. > 1 is an enhancement due to a possible 
increase in the effective degrees of freedom at high den- 
sity m. 

Using eqs. ([5]) and ([S]) one can find the minimal mix- 
ing angle consistent with > 0.25. The minimal 
allowed value of the mixing angle corresponds to the 
maximal allowed value of the ratio of these two tem- 
peratures, Tcoro/Tly_sphere- The neutrinosphcre temper- 
ature T^_sphere = 2-5 McV [13, [H H!] is deter- 
mined by the heat exchange due to neutrino cooling 
near the surface of last scattering. This temperature 
has very little dependence on the nuclear equation of 
state and is almost entirely determined by the conditions 
around the neutrinosphere, i.e. at density of the order of 
10" - 10^2 g/cm^ 24, 25, 26]. In contrast, the core tem- 
perature depends on the nuclear equation of state at den- 
sities p ^ lO^^'g/cm'^ and can vary dramatically, depend- 
ing on the assumptions about nuclear matter [23, HI]. 
With a few exceptions, the models listed in Tables of 
Ref. _26| predict the core temperatures Tcore < 100 MeV, 
and the majority of these models show the temperatures 
in the range Tcore « (20 - 70) MeV We wiU adopt 
the latter as the allowed range of the core temperatures. 

Since there is little or no correlation between the core 
temperature and the neutrinosphere temperature, we al- 
low the ratio Tcore/Ty-sphere Vary between the lowest 
and the highest value for Ty_sphere = 2 — 5 MeV and 
Tcore ~ (20 — 70) MeV. The large uncertainty in the ra- 
tio of temperatures is the reason for the broad allowed 
range of masses and mixing angles. The corresponding 



contours are shown in Fig. 1 for different values of the 
time delay. The longer time delays correspond to lower 
mass for the same mixing angle. However, since the time 
available for the kick is shorter, the minimal mixing angle 
increases for longer delays (see Fig. 1). 

The neutrino-driven kicks have several predictions, in 
addition to time delays, that can be tested using the as- 
tronomical observations. The neutrino kick mechanism 
does not predict a correlation between the magnitude of 
the surface magnetic fields and the pulsar velocity (the 
"S — u" correlation) Q . The kick velocity is determined 
by the magnetic field inside the hot neutron star during 
the first seconds after the supernova collapse. In con- 
trast, the astronomical observations can be used to infer 
the surface magnetic fields of pulsars some millions of 
years later. The relation between the two is highly non- 
trivial because of the complex evolution the magnetic 
field undergoes in a cooling neutron star. However, the 
correlation of the direction of the spin axis and the direc- 
tion of the pulsar velocity is a generic prediction of this 
mechanism Such a correlation is confirmed by 

recent observations [1^ . In the event of a nearby super- 
nova, the neutrino kick can produce gravity waves that 
could be detected by LIGO and LISA [11,13^. Finally, 
the neutrino-driven kicks can increase the energy of the 
supernova explosion because they enhance the convec- 
tion in front of the movin g ne utron star and increase the 
energy of the shock wave [3l|, and also because they de- 
posit entropy ahead of the shock [32] . The increase of 
convection in front of the moving neutron star can pro- 
duce asymmetric jets with the stronger jet pointing in 
the direction of the pulsar motion, in contrast with what 
one could expect from other kick mechanisms f3l|. 
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